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COMMUNICATION DEVICE, COMMUNICATION METHOD, PROGRAM, AND 

STORAGE MEDIUM 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to a communication device, 
a communication method, a program, and a storage medium. 
More particularly, the present invention relates to a 
communication device, a communication method, a program, and 
a storage medium storing a program, which allow multiple 
access in UWB (ultra wideband) communication. 

2 . Description of the Related Art 

UWB communication is also called impulse radio 
communication and is expected to find various applications, 
as discussed for example, in "Impulse radio: How it works" 
by M. Z. Win and R. A. Scholtz (IEEE Communication Letters, 
Vol. 2, pp. 36-38, Feb, 1998 (hereinafter, this literature 
will be denoted as "Ref. 1")), "UWB waveforms ad coding for 
communications and radar" by L Fullerton (Proc. IEEE 
Telesystems Conf., pp. 139-141, Mar. 26-27, 1991 
(hereinafter, this literature will be denoted as "Ref. 2")), 
and "Answers to questions posed by Bob Lucky, Chairman of 
the FCC's Technical Advisory Committee (e-mail 
correspondence, http: //ultra. usc.edu. /ulab, June 29, 1999 
(hereinafter, this literature will be denoted as "Ref. 3")). 
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In UWB communication, data is transmitted using a very 
short pulse at low power without using a carrier. To 
transmit such a very short pulse , a wide bandwidth on the 
order of several GHz is used, as described, for example, in 
Ref. 1 or 2. 

The UWB communication technology has many advantages 
rt that stem from its ultra-wideband nature, as described in, 

a 

,g for example, "On the robustness of ultra-wideband width 

JjJ signals in dense multipath environments" by M. Z. Win and R. 

h£ a. Scholtz (IEEE Communication Letters, Vol. 2, pp. 51-53, 

jL : Feb., 1998 (hereinafter, this literature will be denoted as 

]5f "Ref. 4")), and "On the analysis of UWB communication 

i & 

!i* channels" by J. M. Cramer, R. A. Scholtz, and M. Z. Win 

irs :; 

^ (Proc. IEEE Military Communications Conf., pp. 1191-1195, 

Oct. 31 to Nov. 3, 19 99 (hereinafter, denoted as "Ref. 5")). 
More specifically, it can penetrate walls, experience 
significantly less fading, offer extremely fine time- 
resolution, and deliver large processing gains. 

In particular, the absence of a carrier signal in the 
UWB communication obviates the need for radio frequency (RF) 
or intermediate frequency (IF) circuits. This allows 
reductions in device size and power consumption. 

These characteristics of the UWB communication 
technology give it substantial advantages over conventional 
narrowband, wideband and infrared wireless communication 
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systems . 

Potential applications of the UWB communication include 
wireless local area networks (LAN), information distribution 
systems for use in a medical application or the like, 
outdoor or indoor wireless multiple-access communication 
systems, ranging devices, and other many applications, as 
described in, for example, "Catching the wave; Breakthroughs 
4f in wireless technology" by J. S. Gage (MD Computing: The 
J Leading Edge in Medical and Healthcare Informatics, Vol. 16, 
yJ March/April 1999 (hereinafter, denoted as "Ref. 6")), and 
is "Ultra-wideband width time-hopping spread-spectrum impulse 

yj radio for wireless multiple-access communications" by M. Z. 

W 

Ijfj Win and R. A. Scholtz (IEEE Trans. Communications, Vol. 48, 

.sat. 

}u PP- 679-691, Apr., 2000 ( hereinafter , denoted as "Ref. 7")). 
There is a great need for applying the UWB 
communication technology to multiple-access (multi-user) 
communication with a plurality of communication terminals. 
The idea of applying UWB to multiple-access communication is 
discussed in, for example, "Multiple access with time- 
hopping impulse modulation" by R. A. Scholtz (Proc. IEEE 
Military Communications Conf. (Boston, U.S.), pp. 447 - 450, 
Oct. 11-14, 1993 (hereinafter, denoted as "Ref. 8")), and 
"Multiple-access performance limits with time hopping and 
pulse position modulation" by F. Ramirez -Mireles and R. A. 
Scholtz (Proc. IEEE Military Communications Conf., pp. 529- 



533, Oct. 18-21, 1998 (hereinafter, denoted as "Ref. 9")). 

UWB multiple-access communication systems proposed in 
Refs. 7 to 9 are conceptually similar to asynchronous code- 
division multiple-access (CDMA) systems in certain aspects. 
Transmitted signals of each user (signals transmitted from 
each communication terminal) share a common spectrum, and 
the set of signaling waveforms across all users are not 
necessarily orthogonal. 

The multiple-access capability of a UWB system is 
primarily determined by the processing gain that is given by 
the pulse bandwidth to the pulse (or symbol) repetition 
frequency (RPF) ratio. In the UWB communication, its very 
wide frequency bandwidth results in a high processing gain. 
The UWB communication technology is very different from the 
CDMA communication technology in that the high processing 
gain of the UWB communication allows it to be advantageously 
employed in applications in which power is limited. 

In the literatures Ref. 1 and Ref s . 7 to 9, it is 
proposed to realize multiple-access UWB communication by 
time-hopping the timing of transmitting pulse-position- 
modulated signals of respective users on the basis of a 
time-hopping sequence indicated by a pseudorandom number. 

Herein, the "pulse position modulation" refers to 
modulation of pulse positions depending on data to be 
transmitted. For example, a unit period with a 



predetermined length is assigned to one symbol, and each 
period is divided into M intervals. A pulse is placed in 
one of the M intervals, at a position corresponding to data 
to be transmitted. The "time hopping" refers to an 
operation of, rather than transmitting symbols (unit 
periods) at fixed intervals, randomly shifting the 
transmission timing from the fixed intervals. 

In the proposed UWB multiple-access communication 
techniques, each receiving device includes a single-user 
demodulator for detecting the pulse position, in each unit 
period, of a pulse-position-modulated signal. 

The single-user demodulator includes a filter bank 
including filters for matching the received signal with 
pulses at possible positions in unit periods, and a detector 
for detecting a filter that provides a highest output among 
the filters of the filter bank, and employing, as 
demodulated data, a value corresponding to a pulse position 
that is matched with the received signal by the filter that 
is determined to provide the highest output. 

Ideally, in the multiple-access communication, there is 
no overlap among pulses of signals transmitted from users. 
Under such an ideal condition, the single-user demodulator 
can function as an ideal demodulator for demodulating a 
pulse-position-modulated signal in an environment including 
additive white Gaussian noise. In this case, the single- 
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user demodulator may be constructed in a similar fashion to 
a single-user matched filter for detection in CDMA systems, 
as described, for example, in "Multiuser detection for CDMA 
systems" by A. Duel-Hallen, J. Holtzman and Z. Zvonar (IEEE 
Personal Communications Magazine, Vol. 2, pp. 46-58, Apr., 
j^, 1995 (hereinafter, denoted as "Ref. 10")), "Multi-user 
p detection for DS-CDMA communications", by S. Moshavi (IEEE 
rt Communications Magazine, Vol. 34, pp. 124-136, Oct., 1996 

; rv. 

'rl (hereinafter, denoted as "Ref. 11")), and "Multiuser 

z ~ Detection" by S. Verdu (Cambridge University Press, 1998 

H: (hereinafter, denoted as "Ref. 12")). 

KM 

l r£ However, in practical multiple-access communication 

systems, interference occurs among signals transmitted from 
users (hereinafter, such interference will be referred to as 
multiple-access interference), and multiple-access 
interference (MAI) is generally not Gaussian. The 
conventional single-user modulator is no longer optimal in 
such an environment in which MAI can occur. 

MAI increases with the number of multiple-access users, 
and the increase in MAI causes various adverse effects on 
the single-user demodulator. When data is transmitted at a 
high rate, the resultant high pulse repetition frequency 
causes a reduction in the processing gain and thus an 
increase in MAI. 

In view of the above, it is an object of the present 



invention to provide a high-performance demodulator for use 
in UWB multiple-access communication. 

SUMMARY OF THE INVENTION 

According to an aspect of the present invention , there 
is provided a communication device comprising correlation 
calculation means for calculating the correlation between a 
received signal corresponding to signals transmitted from a 
plurality of communication terminals by means of UWB (ultra 
wideband) communication and pulses at possible positions in 
a signal transmitted from each communication terminal; and 
demodulation means for demodulating data transmitted from 
the respective communication terminals on the basis of the 
correlation, taking into account interference among the 
signals transmitted from the communication terminals. 

According to another aspect of the present invention , 
there is provided a communication method comprising the 
steps of calculating the correlation between a received 
signal corresponding to signals transmitted from a plurality 
of communication terminals by means of UWB (ultra wideband) 
communication and pulses at possible positions in a signal 
transmitted from each communication terminal; and 
demodulating data transmitted from the respective 
communication terminals on the basis of the correlation, 
taking into account interference among the signals 



transmitted from the communication terminals. 

According to still another aspect of the present 
invention, there is provided a program comprising the steps 
of calculating the correlation between a received signal 
corresponding to signals transmitted from a plurality of 
communication terminals by means of UWB (ultra wideband) 
communication and pulses at possible positions in a signal 
transmitted from each communication terminal; and 
demodulating data transmitted from the respective 
communication terminals on the basis of the correlation, 
taking into account interference among the signals 
transmitted from the communication terminals. 

According to still another aspect of the present 
invention, there is provided a storage medium storing 
program comprising the steps of calculating the correlation 
between a received signal corresponding to signals 
transmitted from a plurality of communication terminals by 
means of UWB (ultra wideband) communication and pulses at 
possible positions in a signal transmitted from each 
communication terminal; and demodulating data transmitted 
from the respective communication terminals on the basis of 
the correlation, taking into account interference among the 
signals transmitted from the communication terminals. 

That is, in the communication device, the communication 
method, and the program according to the present invention, 



the correlation is calculated between the received signal 
corresponding to the signals transmitted from a plurality of 
communication terminals by means of UWB (Ultra WideBand) 
communication technology and pulses at possible positions in 
the transmitted signals from the respective communication 
terminals , and the received signal is demodulated into 
original data issued by the communication terminals on the 
basis of the calculated correlation taking into account 
interference among the transmitted signals from the 
respective communication terminals. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram illustrating an example of a 
construction of a UWB communication system according to an 
embodiment of the present invention; 

Fig. 2 is a block diagram illustrating an example of a 
construction of a communication device; 

Fig. 3 is a block diagram illustrating an example of a 
construction of a UWB demodulator; 

Fig. 4 a diagram illustrating a signal transmitted by 
means of UWB communication technology; 

Fig. 5 is a block diagram illustrating an example of a 
construction of a correlator; 

Fig. 6 is a block diagram illustrating an example of a 
construction of a correlator bank; 
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Fig. 7 is a block diagram illustrating another example 
of a construction of a correlator; 

Fig. 8 is a block diagram illustrating an example of a 
construction of a single-user demodulator; 

Fig. 9 is a block diagram illustrating an example of a 
construction of a symbol detector functioning as an optimum 
demodulator; 

Fig. 10 is a flow chart showing the process performed 
by the UWB demodulator functioning as the optimum 
demodulator; 

Fig. 11 is a block diagram illustrating a first example 
of a construction of a symbol detector functioning as a 
suboptimum demodulator; 

Fig. 12 is a flow chart showing the process performed 
by the UWB demodulator functioning as the suboptimum 
demodulator shown in Fig. 11; 

Fig. 13 is a block diagram illustrating a second 
example of a construction of a symbol detector functioning 
as a suboptimum demodulator; 

Fig. 14 is a flow chart showing the process performed 
by the UWB demodulator functioning as the suboptimum 
demodulator shown in Fig. 14; 

Fig. 15 is a graph illustrating a simulation result; 

Fig. 16 is a graph illustrating a simulation result; 

and 



Fig. 17 is a block diagram illustrating an example of a 
construction of a computer according to an embodiment of the 
present invention . 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig. 1 illustrates an embodiment of a UWB communication 
system according to the present invention (herein, the term 
"system" refers to a logical set of apparatuses, the 
physical locations of which may or may not be within a 
single physical enclosure). 

In the embodiment shown in Fig. 1, the UWB 
communication system includes K+1 communication devices l lf 
1 2/ ..., l k+1 . Each communication device l k modulates data to 
be transmitted by means of, for example, pulse position 
modulation and transmits the modulated signal as a radio 
transmission signal using a UWB communication technique 
through a transmission path 2. Each communication device l k 
also receives signals transmitted from the other 
communication devices l k through the transmission path 2 and 
demodulates the received signals into original data. 
Hereinafter, a user of the communication device l k is denoted 
as a user #k. 

In the following discussion, among K+1 communication 
devices 1 1# 1 2 ,..., l k+1/ a communication device l k+1 is taken 
as an example, and the operation thereof is described. The 
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communication device l k+1 is capable of communicating with 
any of the other communication devices l l# 1 2/ ..., l k of K 
users #1, #2,..., tK by means of multiple access 
communication/ in which a received signal corresponding to a 
signal transmitted from each user is demodulated into 
original data issued by the user. 

H' As with the communication device l k+1 , each of the other 

Q 

.sec, 

W communication devices l lf 1 2 ,..., l k is capable of 

O communicating with the other K communication devices by 

\M means of multiple access communication, although the 

communication device l k+1 is taken as an example in the 
yj following discussion. 

.RS5, 

|jf| The UWB communication system such as that shown in Fig. 

jij 1 may be applied to, for example , a wireless LAN system or 
the like. 

Fig. 2 illustrates an example of a construction of the 
communication device l k+1 shown in Fig. 1. Each of the other 
communication devices 1 1# 1 2 ,..., l k may be constructed in a 
similar manner to the communication device l k+1 . 

The communication device l k+1 includes a UWB modulator 11 , 
an antenna 12, and a UWB demodulator 13. 

Data to be transmitted is supplied to the UWB modulator 
11. The UWB modulator 11 modulates the supplied data by 
means of pulse position modulation in which positions of 
pulses (impulses) are determined in accordance with the 



supplied data. The resultant modulated signal is supplied 
to the antenna 12 while performing time hopping. 

Upon receiving the modulated signal from the UWB 
modulator 11, the antenna 12 transmits the modulated signal 
in the form of a radio wave. The antenna 12 also receives a 
radio wave from other communication devices 1 and supplies 
the received signal to the UWB demodulator 13. 

The UWB demodulator 13 demodulates the received signal 
supplied from the antenna 12 in a manner described later and 
outputs the resultant demodulated data. 

Fig. 3 shows an example of a construction of the UWB 
demodulator 13 shown in Fig. 2. 

The UWB demodulator 13 includes a correlator 21 and a 
symbol detector 22. 

The received signal is supplied from the antenna 12 to 
the correlator 21. The correlator 21 calculates the 
correlation of the received signal with respect to possible 
pulse positions for signals transmitted from the respective 
communication terminals l x to l k . The resultant correlation 
is supplied to the symbol detector 22. On the basis of the 
correlation supplied from the correlator 21, and taking into 
account interference among the signals transmitted from the 
communication terminals l x to l k , the symbol detector 22 
demodulates the received signal into original data issued by 
the communication devices l x to l k . 



The process performed by the UWB demodulator 13 shown 
in Fig. 3 is described in further detail below on the basis 
of a model of a multi-user UWB communication system. 

A multiple-access model in an additive white Gaussian 
noise environment has been proposed, for example, in the 
literatures Ref s . 1 and 7 cited above. The multi-user UWB 
communication system can be modeled in part on the basis of 
that given in Ref. 1 or 7 . However, the primary difference 
is the introduction of an equivalent vector model as 
described in "Sequence detection and equalization for pulse- 
position modulation" by J. R. Barry (Proc. IEEE Intl. Conf. 
Communications, pp. 1561-1565, May 1-5, 1994 (hereinafter 
denoted as "Ref. 13")). 

Note that the equivalent vector model can be applied 
not only to the pulse position modulation but also to 
orthogonal modulation schemes. That is, although it is 
assumed that pulse position modulation is employed in the 
present embodiment, the process described below may also be 
applied to orthogonal modulation. 

In the multiple-access model in an environment 
including additive white Gaussian noise, a received signal 
r(t) from user #K + 1 (communication device l k+1 ) can be 
expressed using transmitted signals s <k) (t) from user #k 
(communication devices l k ), as described below: 
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r(r) = JV* } (0 + «(0 (1) 

In this equation (1)/ t denotes time, and n(t) denotes 
zero-mean white Gaussian noise with a power spectral density 
of N Q /2, where N Q is the mean power of white Gaussian noise. 
M' The transmitted signal s (k) (t) of user #k in equation (1) 

can be expressed as follows: 



-CSC" 



□ 



rt 



N M 



^(O-V^FS 2>r m) « m ('-"^-7f>) (2) 



=1 m=l 



wherein E s (k) denotes symbol energy of the transmitted signal 
from user #k. Herein, it is assumed that E s (1) = E s (2) =,..., 
E s <k) = E s/ although the above assumption is not necessarily 
needed, and the invention may also be applied to non-uniform 
s ignal energies . 

In equation (2), N denotes the number of symbols 
transmitted by each user #k (the number of symbols detected 
at a time in detection of a symbol sequence), and T s denotes 
a symbol transmission repetition period for each user #k. 
Thus, 1/T S represents the symbol rate, which is generally 
referred to as the pulse repetition frequency (RPF) in the 
art of UWB communication. 

In equation (2), M denotes the unit symbol size. In 
the case of the pulse position modulation, M indicates the 
number of intervals included in each unit period with a 
predetermined length. Positions corresponding to respective 



M intervals included in each unit period are called chips 
(strictly speaking, M intervals are called chips). In the 
case where each symbol includes M chips and a pulse is 
placed only in one of the M chips depending on data to be 
transmitted, the number of bits transmitted per symbol is 
given by log 2 M. 

In equation (2), b n (k ' m) becomes 1 when a pulse is placed 
at an mth position of an nth symbol of the transmitted 
signal of user #k, and becomes 0 when no pulse is placed at 
that position. Thus, a vector b n k including elements b n (k,m) , 
such as [b n (k,1) , b n (k ' 2) ,..., b n (k ' M) ] T (superscript T represents 
transpose), represents an nth symbol of the transmitted 
signal from the user #k. 

Herein, it is assumed that only one of M elements b n (k ' m) 
in the vector b n <k) is 1, while the remaining M - 1 elements 
are 0. That is, a pulse is placed only at one of M 
positions in each unit period, depending on data to be 
transmitted. More specifically, when j (j=0, 1,...,M- 
1) is data to be transmitted, only (j+l)th element b n <k ' j+1) is 
set to 1, and the other elements are set to 0. Note that, 
alternatively, a plurality of pulses may be placed in each 
unit period depending on data to be transmitted. 

In equation (2), a (m) (t) represents the waveform of a 
pulse at a position m. The pulse a (m) (t) has unit signal 
energy and satisfies the following equation. 
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Sy m \t)a*Xt)dt = 6 m , m . (3) 

where 6 mm is the kronecker delta function. 

More specifically, the pulse a (m) (t) is represented by 

the following equation: 

* {m \t) = w rec (t-mT p ) (4) 

As can be seen from equation (4), the possible 
positions at which pulses a <m) (t) can be placed in each unit 
period are spaced from each other by a time interval T p . In 
other words, the time T p represents time intervals at which 
each unit period is divided into chips (hereinafter, T p will 
be referred to as a chip interval ) . 

The signal w rec (t) that defines the pulse a (m) (t) has a 
very short pulse waveform with a width equal to or smaller 
than 1 ns. For example, as described in Ref. 7, the signal 
w rec (t) can be given by the following equation: 

"«(0 = Jf [1 - 4^r(— ) 2 ]exp[-2 J r(— ) 2 ] ( 5 ) 

' 3 T uwb T uwb 

wherein x uwb is a parameter that determines the pulse width, 
and time t is in units of nano-seconds (ns). For example, 
when x uwb = 0.7 531 ns, the with of the pulse w rec (t) becomes 
about 1.5 ns, as is discussed, for example, in a paper 
entitled "Performance of ultra wideband time-shift-modulated 
signals in the indoor wireless impulse radio channel" by F. 
Ramirez -Mireles , M. Z. Win, R. A. Scholtz and M. A. Barnes 



(Proc. Asilomar Conf . on Signal , Systems and Computers, pp. 
192-196, Nov. 2-5, 1997 (hereinafter denoted as "Ref. 14")) 
and a paper entitled "Performance of Equicorrelated ultra- 
wideband pulse-position-modulated signals in the indoor 
wireless impulse radio channel" by F. Ramirez-Mireles and R 
A. Scholtz (Proc. IEEE Pacific Pim Conf. on Communications, 
Computers and Signal Processings, pp. 640-644, Aug. 20-22, 
1997 (hereinafter denoted as "Ref. 15")). 

Strictly speaking, the pulse width refers to a time 
duration during which the pulse w rec (t) can be regarded as 
existing. 

In equation (5), the factor V(8/3) has been introduced 
to normalize the energy of the pulse w rec (t) such that 



The correlation (autocorrelation) p(t) among pulses 
w rec (t) is given as: 



£.< w ~<0) 2 *-i- 




(6) 



As described, for example, in Ref. 9, the 



autocorrelation p (t) given by equation (5) can be 



approximated as follows: 




(7) 



where a = jt(t/x uwb ) 2 . 
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In equation (2), T n (k) represents a TH (time-hopping) 
time for a symbol of user #k. For example, when symbols are 
transmitted from a user #k to a user #K+1, the timing of 
transmitting symbols is subjected to time hopping in 
accordance with a pseudorandom number. More specif ically, 
an Jith symbol from the user #k is transmitted at a time 
M; shifted from a time nT s by the TH time T n (k) . The time 

O hopping is performed to reduce possible collisions with 

y3 

Q symbols transmitted by other users #k' . 

MS 

U The TH time T n (k) is represented using T <k> and T' n (k) such 

m 

fi that T n <k) = T (k) + T' n (k) . Herein, T (k) can take a value within 
iij the range of 0 ^ T (k) < T s and represents a delay of a symbol 
jjrt transmitted from user #k relative to a symbol delay T (1) for 
j«y user #1 (that is, T (1) = 0). Time T , n (k) represents pseudo 
noise that determines time hopping and can take a value 
within the range of 0 <s T , n (k) < T TH , max . Herein, it is 
required that T TH#max <s T s , and it is also required that the 
minimum difference T s - T TH#max should be smaller than a 
sample-and-hold time of the S/H ( sample-and-hold) circuit 55 
in a correlator 41 m (Fig. 7) that will be described later and 
also than an interval between a pulse placed in a certain 
symbol and a pulse placed at a next symbol. 

Herein, for the purpose of simplicity, it is assumed 
that T TH#max = T s and it is also assumed that the time-hopping 
sequence {T n (k) } is modeled as a sequence of independent and 
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identically distributed continuous random variables each 
uniformly distributed over the range from (equal to or 
greater than) 0 to (smaller than) T THmax . In practice, 
however, the time-hopping sequence {T <k) } are discretized 
over chip intervals • 

As can be seen from equation (2), a signal s <k) (t) is 
transmitted from a user #k as shown in Fig. 4A. 

In Fig. 4A, a symbol #n is transmitted from the user #k 
when a period of time equal to T s + T n (k) has elapsed since 
the transmission of a previous symbol #n-l. 

In the symbol #n, as shown in Fig. 4B, a pulse 
represented by V(E s <k) )w roc (t) is placed at a position m of M 
positions spaced from each other by T p/ wherein the position 
m depends on the data being transmitted. 

Referring again to Fig. 3, the correlator 21 calculates 
the correlation of the received signal r(t) represented by 
equation (1) with respect to possible pulse positions for 
signals transmitted from the respective users #1 to #K. 

For a certain user #k* of the plurality of users #1 to 
#K, a pulse placed at an m'th position in an n'th symbol 
(such a pulse is referred to as a local pulse) of a signal 
transmitted from the user #k' can be represented as a (m ' >(t - 
n'T s - T n . (k ) ). Thus, for the user #k' , the correlator 21 
calculates the correlation y n ,< k, ' m '> i n accordance with the 
following equation: 



=f_j{t)" {m \t-ri*T s -T^)dt (8) 

For the user #k' , the correlation y n ,< k '' m '> represented by 
equation (8) is calculated for all positions m* . Herein, we 
introduce an M x 1 vector y n . (k ' > whose elements are given by 
correlations y n (k "' l) # y n . <k ' ,2) , • . . , y n . (k '' M> calculated for the 
respective positions nT = 1, 2,..., M. That is, the M x 1 
vector y n . <k ) can be represented as follows: 

& ) f (9) 

Because symbol transmission between a certain user and 
another user is performed asynchronously, the correlations 
yn' (1) / y n <2) / • • • /y n (k J associated with the respective users 
outputted from the correlator 21 are not necessarily 
deterministic with respect to the time phase, and the 
correlations should be treated as stochastic variables. 
Therefore, hereinafter, values (such as y n ,< k, ' m *> or y n . <k ) ) 
outputted from the correlator 21 will be referred to as 
correlation statistics . 

The correlation statistic y n ,< k '' m '> i n equation (8) can be 
determined by substituting the received signal r(t) given by 
equation (1) into equation (8). Thus, the correlation 
statistic y n .< k '' m, > i s given as: 

= p »'- k '»' T Ab + V n '' k ' m ' (10) 
where p n, ' k '' m ' is the correlation vector between the local 
pulse a(m')(t - n'T s - T n . (k ) ) and the signals s (1) (t) to s k (t) 
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transmitted from the users #1 to #K, and the correlation 
vector p n *' k '' m ' can be represented in the form of an NKM x 1 
vector as shown below: 

P = [Pi >~>Pn J (11) 

The correlation vector p n, ' k, ' m, in equation (11) includes 
N elements in the form of vectors p n n, ' k '' m ' . Each vector 
p^n',k-,m- ^ s a correlation vector between a local pulse a m '(t - 
n'T s - T n . (k ) ) and an nth symbol in signals s (1) (t) to s (k) (t) 
transmitted from users #1 to #K, and is represented in the 
form of a KM x 1 vector as shown below: 

Pn =[Pn,l >->Pn,K ] ( i2 ) 

The correlation vector p n n ' k ' m in equation (12) includes 
N elements in the form of vectors p n n ' k /Itl • Each vector 
p n n '' k, ' m ' is a correlation vector between a local pulse a (m,) (t - 
n*T s - T n . (k > ) and an nth symbol of a signal s (k) (t) 
transmitted from a user #k, and is represented in the form 
of an M x 1 vector as shown below: 

Pn,K YPnJkA >~->f J n t k t M J V ± J ) 

Each element p n , k , m n,,k '' m ' of the correlation vector p n , k n '' k '' m ' 
in equation (13) represents the correlation between a local 
pulse a (m ) (t - n*T s - T n , (k,) ) and a pulse at an mth position of 
an nth symbol of a signal s(k)(t) transmitted from a user #k, 
and can be calculated in accordance with equation (14) shown 
below, wherein the correlation p n , k , m n ' k ' m ' in this equation 
(14) can be calculated in accordance with equations (4) to 



- 23 - 



'-US 



(7) described above. 

In equation (10), the matrix A is a KMN x KMN diagonal 
matrix representing the signal energies for all symbols and 
all users, and can be represented as follows: 

A-diag^,..^) (15) 
where each element of the matrix A is a KM x 1 vector, and 
diag(A n ) denotes a diagonal matrix whose diagonal elements 
are given by vectors A„. 

Each vector A n that is an element of the matrix A in 
equation (15) can be represented by equation (16) shown 
below. 

\-[^ T ,...,K Kf f (16) 
In equation (16), each element A^ of the vector A„ is a 
M x 1 vector having identical elements E s k as shown in 
equation ( 17 ) . 

=[E?\E?\...,EPf (17) 

In equation (10), the vector b is an NKM x 1 vector 

given by equation (18) representing possible pulse positions 

for a symbol of each user, that is, possible symbol values. 
b f (18) 

Each element b n of the vector b in equation (18) is a KM 
x 1 vector representing possible values of an nth symbol 
transmitted from each user, and is given by equation (19) 
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shown below. 

K=[bf,..,b^f (19) 

Each element b n k of the vector b n in equation (19) is an 
M x 1 vector representing possible values of an nth symbol 
transmitted from a user k, and is given by equation (20) 

shown below. 

bj» ~[bf^ t .„ t tf^f (20) 

In equation (10), T) n, ' k,,m ' is a Gaussian random variable 
with a mean value of 0 and a variance of N G /2. 

If the collection of correlation statistics y n (k > for 
all symbols and all users, which are represented by equation 
(9) and outputted from the correlator 21 , is denoted by a 
vector y, then equation (10) can be rewritten using this 
correlation statistic y as follows: 

y = RAb + 77 (21) 
where the vector b is defined in equations (18) to (20). 

The vector y in equation (21) is a KMN x 1 vector 
including N vectors given by the following equation: 

y = [yl^,ylf (22) 

In equation (22), each element of the vector y, that is, 
the vector y n . , represents a correlation statistic of an iith 
symbol for all users, and the vector y n . is a KM x 1 vector 
represented by the following equation: 

y n > = [yf,~,yP T f (23) 

In equation (23), the vector y n . includes elements in 
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the form of vectors y n (k<) defined in equation (9). 

In equation (21), R is an NKM x NKM correlation matrix, 
which is defined as follows : 

R-[p lT 9 ... 9 p» T f (24) 

In equation (24) , each element p n ' of the correlation 
matrix R is a KM x NKM matrix representing the correlation 
between a local pulse and an n'th symbol for all users, and 
is given by the following equation: 

p =[ p , , p ] (25) 

In equation (25) , each element p n ' k '' of the matrix p n ' is 
an M x NKM matrix representing a local pulse and an n'th 
symbol of a signal transmitted from a user #k' , and is given 
by the following equation: 

P =[ P 9 • • • t P ] (26) 

In equation (26) , each element p n '' k, ' m, Q f the matrix p n ' k *' 
is defined in equation (11). 

r| in equation (21) corresponds to y ^ n '' k, ' m, i n equation 
(10) and defined by a vector (noise vector) given by the 
following equation : 

rj = [rtl , .,rj T N ] T (27) 
where the vector rj n . is defined as follows: 

•fc- = ir,f ,...,r,^ T V (28) 
where the vector Ti n . (k > is defined as follows: 



in /•••/*/ j 



(29) 
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Note that elements of r] n k ' in equation (29) are given by 
^n-.k'..- in equation (10). 

The correlator 21 shown in Fig. 3 calculates the 
correlation statistic y in equation (21) (in practice, the 
correlator 21 calculates the correlation statistic in 
accordance with equation (8), which results in the same 
correlation statistic according to equation (21). To this 
end, as shown in Fig. 5, the correlator 21 includes as many 
K correlator banks 31 x to 31 k as there are multi-access users 
and also includes switches 32 x to 32 k disposed at the outputs 
of the respective correlator banks 31 1 to 31 k . 

The received signal r(t) is supplied to each correlator 
bank 31 k .. The correlator bank 31 k . calculates the correlator 
statistic y n <k > in equation (9) using the received signal 
r(t) and the local pulse a (m ) (t - n'T s - T n . <k ) ) in accordance 
with equation (8). 

The correlation statistic y n . (k>) outputted from the 
correlator bank 31 k . is supplied to the switch 32 k .. The 
switch 32 k . is normally in an off-state and is turned on at a 
time t = (n 1 + 1)T S so that the correlation statistic y n (k ) 
outputted from the correlator bank 31 k . is supplied to the 
symbol detector 22 (Fig. 3) disposed at the following stage 
via the switch 32 k . at t = (n 1 + 1)T S . 

Fig. 6 shows an example of a construction of the 
correlator bank 31 k . shown in Fig. 5. 



The correlator bank 31 k . includes as many correlators 41 x 
to 41 m as there are M possible pulse positions within a unit 
period during which a symbol is transmitted. 

The received signal r(t) is supplied to each correlator 
41 ro . . The correlator 41 m . calculates the correlator statistic 
y n . <k,/m ) in equation (8) using the received signal r(t) and 
the local pulse a (m ) (t - n'T s - T n , (k ) ) in accordance with 
equation ( 8 ) . 

Fig. 7 shows an example of a construction of the 
correlator 41 m< shown in Fig. 6. 

A TH (time-hopping) time generator 51 generates a TH 
time T n , (k,) and supplies it to a pulse generator 52. On the 
basis of the TH time T n . (k ) supplied from the TH time 
generator 51 , the pulse generator 52 generates a local pulse 
a (m,) (t - n'T s - T n . (k ) ) and supplies it to a multiplier 53. 

In addition to the local pulse a (m ) (t - n'T s - T n . (k ) ) 
supplied from the pulse generator 52, the received signal 
r(t) is also supplied to the multiplier 53. The multiplier 
53 calculates the product of the received signal r(t) and 
the local pulse a (m ) (t - n f T s - T n . (k ) ) thereby determining 
r(t)a (m ) (t - n'T s - T n . (k ) ) that is to be integrated with 
respect to time as shown in equation (8). The result is 
supplied to an integrator 54. 

The integrator 54 integrates the output of the 
multiplier 53 with respect to time in accordance with 
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equation (8). The resultant integral value to an S/H 
circuit 55. 

The S/H circuit 55 samples and holds the output of the 
integrator 54 at a time indicated by the TH time T n , (k ) . The 
sampled and held value is outputted as the correlation 
statistic y n ,< k, ' m '> represented in equation (8). 

In the case where the UWB demodulator 13 shown in Fig. 
3 is constructed as a single-user demodulator, the 
correlator 21 includes, as shown in Fig. 8, only one 
correlator 61 similar to the correlator bank 31 k . shown in 
Fig. 6. 

In this case, the output of the correlator 21 
represents the correlation statistic y n> < k '' m, > (= [y n .. (k,fl) , 
y n . (k '' 2) , . . . , y n (k ,M) ]) associated with only a particular user 
#k' , and the symbol detector 22 at the following stage 
detects a symbol on the basis of the correlation statistic 
y n . (k ' m ) associated with that user #k* . 

In this case, because only the correlation statistic 
y n . (k ' m ) of the single user #k' is supplied to the symbol 
detector 22, it is impossible to take into account the 
correlation statistics for other users, that is, 
interference with other users, in the symbol detection. 
That is, the symbol detector 22 estimate a value b , n . (k,,m ' ) for 
a symbol from a user #k' without taking into account 
interference with other users, in accordance with a 



determination criterion given by the following equation. 



That is, in the case where the UWB demodulator 13 is 
constructed as the single-user demodulator, the symbol 
detector 2 2 detects a maximum value y n . <k *' m < max >> from the 
correlation statistics y n . (k ' fl) , . y n . (k ' ,2) , . . . , y n ' (k ' M> # and m(max) 
corresponding to the detected maximum value is outputted as 
demodulated data. 

In contrast, in the case where the UWB demodulator 13 
is constructed as an optimum multi-user demodulator (optimum 
detector), the symbol detector 22 may be constructed, for 
example, as shown in Fig. 9. 

In the example shown in Fig. 9, the symbol detector 22 
includes a correlation matrix calculator 71, an energy 
estimator 72, a symbol pattern generator 73, and a maximum- 
likelihood sequence detector 74. 

The correlation matrix calculator 71 calculates the 
correlations among pulses at possible positions for signals 
transmitted from users in accordance with equation (14) 
thereby determining the correlation matrix R including the 
calculated correlations as represented in equation (24), and 
supplies the result to the maximum- likelihood sequence 
detector 74. In the above calculation of the correlation 
matrix R performed by the correlation matrix calculator 71 




(30) 
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in accordance with equation (14), time information such as T p/ 
T s , and T n k , and the number M of intervals in a unit symbol 
period is needed. The information about these (correlation 
matrix calculation information) is assumed to be given to 
the correlation matrix calculator 71. 
M 1 The energy estimator 72 estimates the energies E s <k) of 

the transmitted signals for respective users #k and 
determines the matrix A (energy matrix) including the 
Uf estimated energies E s k as elements as shown in equation (15). 



m 



The resultant energy matrix A is supplied to the maximum- 
likelihood sequence detector 74. In the present embodiment , 



yfl as described above, the energies E s ( E s ( E s ( ' of the 

□ 

Fy signals transmitted from the respective users are assumed to 
have an identical value E s . 

The symbol pattern generator 73 determines possible 
pulse positions for symbols from the respective users. That 
is, the symbol pattern generator 73 determines the vector b 
(pattern vector), which is defined in equation (18) and 
which represents a possible symbol pattern. The resultant 
pattern vector is supplied to the maximum-likelihood 
sequence detector 74. 

To the maximum- likelihood sequence detector 74, the 
correlation statistic (vector) y outputted from the 
correlator 21 is also supplied in addition to the 
correlation matrix R outputted from the correlation matrix 



calculator 71, the energy matrix A outputted from the energy 
estimator 72, and the pattern vector b outputted from the 
symbol pattern generator 73, as described earlier. Using 
information about these, the maximum-likelihood sequence 
detector 74 determines the maximum- likelihood sequence 
thereby determining an estimated symbol sequence (vector) 
b ,(1) , b ,(2) ,..., b ,(k) and outputs the resultant estimated 
symbol sequence as demodulated data. 

More specifically, the maximum- likelihood sequence 

o 

detector 2 2 determines a maximum- likelihood symbol sequence 
(vector) b 1 = [b ,(1) , b ,(2) ,..., b ,(k) ] for the users #1 to #K, 
using the correlation statistic y, the correlation matrix R, 
the energy matrix A, and the pattern vector b, in accordance 
with a determination criterion given by the following 
equation. 

V = arg min||}> - RAbf (31) 

in equation (31), arg min(x) with respect to b 
represents an operation of determining b that results in 
minimization of x. 

In accordance with equation (31), of M KN pattern vectors 
b, a vector that gives a minimum value of the norm ||y - RAb|| 2 
is employed as the moist-likelihood symbol sequence b 1 . 

In the case where the UWB demodulator 13 is constructed 
as the optimum multi-user demodulator as shown in Fig. 9, 



the UWB demodulator 13 performs a process in accordance with 
a procedure shown in the form of a flow chart in Fig. 10. 

First, in step SI, the correlation matrix calculator 71 
determines the correlation matrix R, and the energy 
estimator 72 determines the energy matrix A. Furthermore, 
the symbol pattern generator 73 determines the pattern 
vector b and supplies the result to the maximum-likelihood 
sequence detector 74. 

Thereafter, in step S2, the correlator 21 determines 
the correlation between the received signal r(t) and the 
local pulses. The correlation statistic y obtained as a 
result is supplied to the maximum-likelihood sequence 
detector 74. 

Using the correlation statistic y, the correlation 
matrix R, the energy matrix A, and the pattern vector b, the 
maximum- likelihood sequence detector 74 determines the 
maximum-likelihood sequence b 1 = [b ,(1) , b ,<2) ,..., b ,(k) ] for 
the users #1 to #K in accordance with a determination 
criterion given by equation (31) and outputs the result as 
the demodulated data. 

Thereafter, the process returns to step S2 to perform 
the above-described operation for the next received signal 
r(t). 

Although the UWB demodulator 13 can be constructed in 
the form of an optimum multi-user demodulator (optimum 
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detector) as described above, the UWB demodulator 13 may 
also be constructed in the form of a suboptimum multi-user 
demodulator (suboptimum detector). 

That is, the UWB demodulator 13 may detect optimum 
symbols (sequence) from the respective users by jointly 
demodulating the signals from all users on the basis of all 
correlation statistics y n . (1) , y n » <2) /***/ y n <k> °f all users 
taking into account interference among users, or may detect 
suboptimum symbols (sequence) by removing interference 
components among users (multiple-access interference) from 

the correlation statistics y n . (1) , y n (2) / / yn (k) using the 

correlation statistics y n . (1) , y n <2) f*-*/ Yn (k) of a H users. 

The UWB demodulator 13 may be constructed in the form 
of a suboptimum multi-user demodulator, for example, using a 
decorrelator . 

The decorrelator is analogous to a zero-forcing linear 

detector for intersymbol interference limited channels in 

single-user systems. The decorrelator projects the 

correlation statistic y onto a signal space having no 

multiple-access interference (MAI) by means of a linear 

transformation in accordance with equation (32) shown below 

thereby determining the correlation statistic y in the 

signal space. 

y=R~ l y (32) 

After determining the correlation statistic y' in the 



signal space having no multiple-access interference in the 
above-descried manner, a suboptimum symbol sequence b' n , (1 ' m) , 

b' n . (2 ' m) / b' n . (k#m) for the users #k to #K can be estimated in 

accordance with the determination criterion given by 
equation ( 30 ) . 

Fig. 11 shows an example of a construction of the 
symbol detector 22 in the UWB demodulator 13 constructed as 
a suboptimum multi-user demodulator using a decorrelator • 

As with the correlation matrix calculator 71 shown in 
Fig. 9, a correlation matrix calculator 81 calculates the 
correlation matrix and supplies the result to an inverse 
matrix calculator 82. The inverse matrix calculator 82 
calculates the inverse matrix R" 1 of the correlation matrix R 
supplied from the correlation matrix calculator 81 
(hereinafter, such an inverse matrix will be referred to as 
an inverse correlation matrix) and supplies the result to an 
inverse correlation calculator 83. 

If the determinant of the correlation matrix R is equal 
to 0 or a very small value, which can occur when there is 
noise or a calculation error or when the correlation matrix 
R is not regular, the inverse matrix calculator 82 performs 
exception handling in which the inverse correlation matrix R" 
1 is replaced by a generalized inverse matrix with a maximum 
square norm. 

In addition to the inverse correlation matrix R" 1 
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outputted from the inverse matrix calculator 83, the 
correlation statistic y outputted from the correlator 21 is 
also supplied to the inverse correlation calculator 83. The 
inverse correlation calculator 83 converts the correlation 
statistic y into a correlation statistic y' in a signal 
space including no multiple-access interference, in 
accordance with equation (32), that is, by means of a linear 
transformation using the inverse correlation matrix R" 1 . The 
resultant correlation statistic y f is supplied to a maximum 
value detector 84. 

In accordance with the determination criterion given by 
equation (30), the maximum value detector 84 determines the 
maximum value y ■ n . < k, ' m < max >> Q f the correlation statistics 
y* nl (k,,1) , y , n . (k '' 2) , . . . , y n ,< k '' M > associated with the user #k in 
the correlation statistic y* supplied from the inverse 
correlation calculator 83. The maximum value detector 84 
then determines m(max) at which the maximum value 
y • ni <*-,m(max)) is obtained. Furthermore, on the basis of the 
determined m(max), the maximum value detector 84 estimates a 
suboptimum symbol sequence b' n . (k,) for the user #k and outputs 
it as demodulated data. The maximum value detector 84 also 
determines and outputs demodulated data for other users in a 
similar manner. 

In Fig. 11, the correlation matrix calculator 81, the 
inverse matrix calculator 82, and the inverse correlation 



calculator 83 form the decorrelator . 

Fig. 12 is a flow chart showing the process performed 
by the UWB demodulator 13 constructed in the form of the 
suboptimum multi-user demodulator using the decorrelator 
such as that shown in Fig. 11. 

In the first step Sll, the correlation matrix 
calculator 81 determines the correlation matrix R and 
outputs the result to the inverse matrix calculator 82. 
Thereafter, the process proceeds to step S12. In step S12, 
the inverse matrix calculator 82 determines the inverse 
matrix R" 1 of the correlation matrix R and supplies the 
resultant inverse matrix R" 1 to the inverse correlation 
calculator 83. 

Thereafter , the correlator 21 determines the 
correlation between the received signal r(t) and the local 
pulses and supplies the correlation statistic y obtained as 
a result to the inverse correlation calculator 83. In step 
S13, in accordance with equation (32), the inverse 
correlation calculator 83 determines the correlation 
statistic y f by mapping the correlation statistic y into a 
signal space including no multiple-access interference. The 
result is outputted to the maximum value detector 84 . 

In step S14, the maximum value detector 84 determines 
the maximum value of the correlation statistic y' supplied 
from the inverse correlation calculator 83 for each user in 



- 37 - 

accordance with the determination criterion given by 
equation (30). Furthermore, on the basis of the maximum 
value, the maximum value detector 84 estimates a suboptimum 
symbol sequence b* for each user and outputs it as 
demodulated data. 

Thereafter, the process returns to step S13, and the 
above-described operation is performed again when the 
correlation statistic y associated with the next received 
signal r(t) is supplied from the correlator 21. 

In the case where the UWB demodulator 13 is constructed 
in the form of an optimum demodulator, as described above 
with reference to Figs. 9 and 10, correlation matrix 
information is required when the correlation matrix R is 
calculated, and estimated energy of a signal transmitted 
from each user is required when the energy matrix A is 
calculated. In the case where the UWB demodulator 13 is 
constructed in the form of a suboptimum demodulator using a 
decorrelator as shown in Fig. 11, calculation of the energy 
matrix A is not needed although calculation of the 
correlation matrix R is required. That is, in this case, it 
is not needed to estimate the energy of the transmitted 
signal from each user. 

When the WUB demodulator 13 is constructed using the 
decorrelator, although the correlation statistic y is mapped 
by the decorrelator into the signal space including no 
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multiple-access interference and thus interference 
components are removed, a slight increase in noise occurs. 

Although, in the above example, a decorrelator is used 
to construct the UWB demodulator 13 into the form of a 
suboptimum multi-user demodulator (suboptimum detector), an 
MMSE (Minimum Mean Square Error) detector may also be 
employed . 

The MMSE detector converts the correlation statistic y 
such that the mean square error thereof becomes minimum. 
This allows well-balanced reductions in intersymbol 
interference by means of zero-forcing method and white 
Gaussian noise. For example, the MMSE detector for this 
purpose may be constructed in a similar manner to a single- 
user linear MMSE detector. 

After conversion of the correlation statistic y using 
the linear MMSE detector, a suboptimum symbol sequence 
b' n . (1 ' m) , b f n . (2 ' m) , . . . b' n . <k ' m) for the users #k to #K can be 
estimated in accordance with the determination criterion 
given by equation (30). 

Fig. 13 shows an example of a construction of the 
symbol detector 22 in the UWB demodulator 13 constructed as 
a suboptimum multi-user demodulator using a linear MMSE 
detector. 

As with the correlation matrix calculator 71 shown in 
Fig. 9, a correlation matrix calculator 91 calculates the 
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correlation matrix and supplies the result to a mean square 
error calculator 93. An energy estimator 92 determines the 
energy matrix A, as with the energy estimator 72 shown in 
Fig. 9, and supplies it to the mean square error calculator 
93. 

In addition to the correlation matrix R outputted from 
the correlation matrix calculator 91, the energy matrix A 
outputted from the energy estimator 92 and the correlation 
statistic y outputted from the correlator 21 are also 
supplied to the mean square error calculator 93. The mean 
square error calculator 93 converts the correlation 
statistic y outputted from the correlator 21 into y' with a 
minimized mean square error using a conversion matrix 
expressed using the correlation matrix R and the energy 
matrix A as shown in equation (33) and supplies the result 
y' to the maximum value detector 94. 

y'-[R+^(AAy l r l y (33) 

In accordance with the determination criterion given by 
equation (30), the maximum value detector 94 determines, as 
with the maximum value detector 84 shown in Fig. 11, the 
maximum value y 1 n< < k, ' m < max >> Q f the correlation statistics 

y , n' (lc ' 1) / y , n (k ,2> / f y n .< k, '*> associated with the user #k in 

the correlation statistic y' supplied from the mean square 
error calculator 93. Furthermore, on the basis of a m(max) 
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at which the maximum value y' n , (lt, ' m(max)) is obtained, the 
maximum value detector 94 estimates a suboptimum symbol 
sequence b' n . <k J for the user #k' and outputs it as 
demodulated data. The maximum value detector 94 also 
determines and outputs demodulated data for the other users 
in a similar manner. 

In Fig. 13 , the correlation matrix calculator 91, the 
energy estimator 92, and the mean square error calculator 93 
form the linear MMSE detector. 

Fig. 14 is a flow chart showing the process performed 
by the UWB demodulator 13 constructed in the form of the 
suboptimum multi-user demodulator using the linear MMSE 
detector such as that shown in Fig. 13. 

First, in step S21, the correlation matrix calculator 
91 determines the correlation matrix R, and the energy 
estimator 92 determines the energy matrix A. The resultant 
correlation matrix R and the energy matrix A are supplied to 
the mean square error calculator 93. 

The correlator 21 determines the correlation between 
the received signal r(t) and the local pulses, and the 
correlation statistic y obtained as a result is supplied to 
the mean square error calculator 93. In step S22, the mean 
square calculator 93 converts the correlation statistic y 
into a correlation statistic y 1 in accordance with equation 
(33) and supplies it to the maximum value detector 94. 



In step S23, the maximum value detector 94 determines 
the maximum value of the correlation statistic y' supplied 
from the mean square error calculator 93 for each user in 
accordance with the determination criterion given by 
equation (30). Furthermore, on the basis of the maximum 
value, the maximum value detector 94 estimates a suboptimum 
symbol sequence b' for each user and outputs it as 
demodulated data. 

Thereafter, the process returns to step S22, and the 
above-described operation is performed again when the 
correlation statistic y associated with the next received 
signal r(t) is supplied from the correlator 21. 

Various types of cancellers for canceling interference, 
such as parallel interference cancellers, successive 
interference cancellers, and decision feedback cancellers, 
are described, for example, in the literatures Refs. 10 to 
12 or in "Analysis of a simple successive interference 
cancellation scheme in a DS/CDMA system", by P. Patel and J. 
M. Holtzman (IEEE J. Select Areas Communications, Vol. 12, 
pp. 796-807, June 1994 (hereinafter, denoted as Ref. 16). 
The principle behind their design is the construction of MAI 
and its subtraction from the received signal to generate 
improved symbol decisions on the basis of the correlation 
matrix. In the case where the UWB demodulator 13 is 
constructed in the form of a suboptimum multi-user 
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demodulator ( suboptimum detector), one of such interference 
cancellers cited above or any combination of such 
interference cancellers may be employed. 

That is, in the case were the UWB demodulator 13 is 
constructed in the form of a suboptimum multi-user 
demodulator using an interference canceller, the 
interference canceller 23 may be disposed between the 
correlator 21 and the symbol detector 22 as represented by a 
broken line in Fig. 3. 

In this case, using all correlation statistics y n . (1> , 
Yn' <2) /---' Yn« <k) outputted from the correlator 21, the 
interference canceller removes interference components among 
users (multiple-access interference) from the correlation 
statistics y n . (1) , y n (2) / — -/ Yn (k) associated with the 
respective users, and outputs the result to the symbol 
detector 22. The symbol detector 22 estimates a suboptimum 
symbol (sequence) for each user, in accordance with the 
determination criterion given by equation (30). 

As described above, in the present embodiment, the UWB 
demodulator 13 calculates the correlation between the 
received signal r(t) corresponding to the signals 
transmitted from a plurality of users by means of the UWB 
communication scheme and pulses (local pulses) at possible 
positions of the transmitted signals from the plurality of 
users. On the basis of the obtained correlation statistic y, 
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and taking into account interference among the signals 
transmitted from the plurality of users, the received signal 
r(t) is demodulated into original data issued by the 
plurality of users. Thus, a high-performance UWB 
communication demodulator is achieved according to the 
present embodiment of the invention. 

For example, in the DS-CDMA (Direct Sequence CDMA) 
communication, the processing gain determined by the ratio 
of the chip rate to the data rate is limited so that the 
correlation among signals transmitted from users does not 
become equal to 0. That is, in the CDMA communication, at a 
transmitting device, data is multiplied by a spread code in 
one period thereby spreading the spectrum of the data over a 
wide band. At a receiving device, a received signal is 
multiplied by the same spread code as that employed at the 
transmitting device and the resultant product is integrated 
over one period of the spread code thereby determining the 
correlation with the received signal. Thus, in the CDMA 
communication, the correlation between spread codes assigned 
to different users remains for a long time. 

In contrast, in the UWB communication, the processing 
gain is determined by the bandwidth-to-PRF ratio. Because 
the bandwidth is very wide as described above, the 
processing gain becomes large. Thus, there is a high 
probability that the correlation among signals transmitted 
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from different users becomes 0. That is, in the UWB 
communication, no spread code is used, but data is 
transmitted using very short pulses which cause the 
frequency spectrum to be distributed over an ultra wide band. 
Therefore, in the UWB communication, the correlation in the 
time domain is localized. Thus, in the UWB communication, 
the correlation matrix R becomes a sparse matrix (most 
elements of which are 0). In the correlation matrix R, non- 
zero elements representing correlations indicate that a 
collision occurs between pulses, that is, multiple-access 
interference occurs. The non-zero correlation increases 
with increasing M indicating the number of chip intervals in 
a unit symbol period. A reduction in the processing gain 
and an increase in the number of multiple-access users K 
(user capacity) also result in an increase in the non-zero 
correlation. 

The dimension of the correlation matrix R increases in 
proportion to the number of users K and the number of 
symbols N of a symbol sequence to be estimated. The 
increase in the dimension of the correlation matrix R 
results in a great increase in the complexity of the 
calculation performed by the UWB demodulator 13. 

However, if the correlation matrix R is sparse, the 
correlation matrix R can be divided into lower-dimensional 
matrices in the calculation. This allows a reduction in the 



complexity of the calculation performed by the UWB 
demodulator 13. For example, in the case where the UWB 
demodulator 13 includes the interference canceller 23 as 
shown in Fig. 3, the interference canceller 2 3 is needed to 
perform a canceling process only for users among which 
multiple-access interference can occur on the basis of the 
TH time T n (k) . This allows a reduction in the complexity of 
the process. 

When the condition KMT p < T s is satisfied, symbols #1, 
#2,..., #N in a symbol sequence do not have an overlap, at 
rising-up/falling-down edges, between pulses corresponding 
to symbols. This makes it unnecessary to deal with a 
problem caused by overlapping between pulses, and thus, the 
UWB demodulator 13 is needed to detect not a sequence of 
symbols but a symbol. That is, when the symbol period 
satisfies the condition MT p < T s , any symbol (pulse) of any 
user is not interfered with by nearby symbols, and thus the 
correlation statistics associated with the respective users, 
outputted from the correlator 21, become independent of each 
other. 

In the UWB communication, as described above, unlike 
the CDMA communication in which the correlation among 
signals transmitted from different users does not generally 
become 0, the correlation matrix R becomes sparse, and thus 
a drastic reduction in the complexity of the process can be 
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In equation (35) , each matrix r n , k n '' k '' is an M x M matrix 



-n, k,m 



n ' , k * ,n 



consisting of elements each representing the 



correlation between a pulse a (m) (t - nT s - T n (k) ) located in an 
nth symbol of a user #k in the received signal r(t) and an 
n'th local pulse a (m ) (t - n'T s - T n . (k ) ) of a user #k' , and the 
matrix r n/k n '' k '' can be represented by M x M matrices r n , k , m n ' ,k '' m ' 
as shown below. 
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(36) 



Of the matrices r n n ' given by equation (35), those 
matrices whose n and n' satisfy the condition |n' - n| > 1 
are zero matrices. Therefore, the correlation matrix R 
given by equation (34) can be represented into the form of a 
matrix including N x N matrices as elements, as shown below: 
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In equation (37), matrix O x denotes an x x x zero matrix, 
From equation (37), the correlation matrix R becomes a 
triangular symmetric matrix similar to that described, for 
example, in "Matrix Computations" 3rd ed., by G. H. Golub 



and C. F. Van Loan (The Johns Hopkins University Press , 1996, 
pp. 16-17, (hereafter, denoted as "Ref. 17")). 

Herein, the "triangular symmetric matrix" refers to a 
matrix whose elements are 0 except for diagonal elements and 
adjacent elements on both sides of diagonal elements. When 
the correlation matrix R is expressed using matrices r n n * 
given by equation (35) as elements, the correlation matrix R 
becomes a triangular symmetric matrix. In other words, if 
the correlation matrix R is represented using elements which 
are all scalar, it does not become a triangular symmetric 
matrix. In this regard, the correlation matrix R is not a 
triangular symmetric matrix in the strict sense. 

When k = k f , elements correspond to symbols from the 
same user. Because there is no correlation among nearby 
symbols from the same user, the matrix r n/k n '' k '' (= r n k n ' k ) 
given by equation (36) becomes a zero matrix 0 m . In the case 
where k = k' and n = n', elements correspond to the symbol 
from the same user. In this case, because pulses at 
different possible positions for the same symbol are 
orthogonal to each other, the correlation p n , k , m n , ' k '' m ' ( = 
Pn,k,m n km ) i n equation (14) can be calculated as follows: 



Therefore, when k = k 1 and n = n' , the matrix r n , k , n ' ,k '' 
given by equation (36) becomes an M x M unit matrix I m . 




(38) 
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Herein, I x denotes an x x x unit matrix. 

Herein, let us assume that the number of users, K, is 
equal to 3, and symbol detection (N = 1) is performed for a 
received signal modulated by means of binary pulse position 
modulation (that is, M = 2). 

Let us further assume that the chip interval T p in 
equation (4) is equal to 1.5 ns, the parameter x uwb in 
equation (5) is equal to 0.7531, and the symbol transmission 
interval T s in equation (2) is greater than 10T p . 
Furthermore, as for the TH time T n k in equation (2), it is 
assumed that T^ 1 * = 0, T 1 (2) = T p , and T x (3) = 5T p . Still 
furthermore, the energy of the transmitted signal is 
identical for all three users #1 to #3. 

Under the above conditions, the correlation matrix R 
can be determined as a6 (=1x3x2) x6 matrix as 
follows : 
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(39) 



Herein, symbols transmitted from the users #1 to #3 can 
be determined by means of the maximum-likelihood estimation 
in accordance with equation (31). That is, in accordance 
with equation (31), the estimated symbols (vector) b 1 is 
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given by a 6 x 1 vector as shown below. 

b'=[b< l >\bf? ,bf?f (40) 

In equation (40), each element b' n k of the vector b 1 is 
a 2 (= M) x 1 vector. 

In the case where the correlation matrix R is given by 
equation (39), the problem of detecting symbols from the 
three users #1 to #3 reduces to detecting two symbols as 
described below. 

That is, the correlation matrix R in equation (39) can 
be divided into four matrices : a 4 x 4 matrix on the upper 
left side; a 2 x 4 matrix on the lower left side; a 4 x 2 
matrix on the upper right side; and a 2 x 2 matrix on the 
lower right side, as shown in equation (41). This means 
that although multiple-access interference occurs between 
the user #1 and the user #2, the user #3 has no multiple- 
access interference with the other users. 
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Herein, if the 4x4 matrix on the upper right side in 
equation (41) is denoted by R' and the x x y zero matrix is 
denoted by O xy , then the correlation matrix R in equation 
(41) can be expressed as follows: 
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(42) 



From equation (42), the detection of symbols according 
to equation (31) reduces to detecting symbols from the user 
#1 and the user #2 between which multiple-access 
interference occurs and detecting the symbol from the user 
#3 having no interference with the other users. 
|SJ Thus, the estimated symbols [b^* 1 *, b\ {2) ] from the user 



#1 and the user #2 can be determined on the basis of 
equation (43) corresponding to equation (31), and the 
estimated symbol b 1 ^ 35 from the user #3 can be determined on 
the basis of equation (44) corresponding to equation (31). 



S [M 1)r ,M 2)r f=arg min r [yf,y[ 2f f -R'W ,t^f 



(43) 



^ = argmin||y< 3 )-/A (3) f (44) 

In the UWB communication, because the correlation 
matrix R has a sparse form as described above with reference 
to, for example, equation (39), the correlation matrix R can 
be divided (reduced) as shown in equation (42), and thus the 
calculation cost can be drastically reduced compared with 
the cost needed when the correlation matrix R is not (cannot 
be) divided. 

In the case where in all N symbols from all users #1 to 
#K, pulses at M possible positions are orthogonal to each 
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other, that is, in the case where there is no overlap among 
pulses, the correlation matrix R becomes an NKM x NKM unit 
matrix I n , k , m . When the correlation matrix R becomes a unit 
matrix as in the above-described case, the process performed 
by the multi-user UWB demodulator 13 becomes equivalent to 
the process performed by K single-user demodulators. 

Figs. 15 and 16 shows the simulated value of the SER 
(Symbol Error Rate) as a function of E b /N Q for the UWB 
demodulator 13 constructed in the form of an optimum multi- 
user demodulator (optimum detector). 

The simulation shown in Figs. 15 and 16 was performed 
by means of Monte-Carlo method on the assumption that M = 2, 
T p = 1.5 ns, x uwb = 0.7531, and the signal energy from any 
user has an identical value E s . 

In Fig. 15, the number of users, K, is assumed to be 2 , 
and the TH times T 1 il) and T x {2) are assumed to be 0 and T Pf 
respectively. Furthermore, in Fig. 15, the symbol 
transmission period T s is set so as to satisfy the condition 
T s > KMT p (= 4T p ) and thus so that N = 1, that is, a symbol 
is detected. The bit energy E b is calculated according to E b 
= E s /log 2 M. 

In Fig. 15 (and also in Fig. 16), a curve labeled "ML" 
indicates the SER for the optimum multi-user demodulator, 
and a curve labeled "CONV" indicates the SER for the single- 
user demodulator. In the single-user demodulator, as can be 
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seen from Fig. 15 , no significant improvement in the SER is 
obtained if the E b /N 0 is increased , and an error floor due to 
multiple-access interference occurs. That is, in the 
single-user demodulator/ because multiple-access 
interference is not taken into account, a pulse placed in 
the first symbol to transmit data #d from a user #1 and a 
pulse placed in the first symbol to transmit the same data 
#d from a user #2 are both detected by the demodulator of 
the user #1 as a pulse placed in the first symbol from the 
user #1, even in an ideal environment including no noise. 
That is, because M = 2 in this example, the unit symbol 
interval is divided into two chip intervals, and a pulse is 
placed only in one of the two chip intervals. Herein, the 
TH times for the user #1 and user #2 are respectively 0 and 
T p , and thus the difference between them is equal to the chip 
interval T p . As a result, the demodulator of the user #1 
detects pulses at both positions in the first symbol. In 
fact, the simulation indicates that the correlation 
statistic has the same value for both positions. The 
detection by the demodulator of the user #2 becomes similar 
for the second symbol transmitted from the users #1 and #2. 

In contrast, in the optimum multi-user demodulator, 
because the multiple-access interference is taken into 
account, the SER is improved as the E b /N 0 increases. 

In Fig. 15, a curve labeled "SINGLE-USER" indicates the 
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lower limit of the SER of the single-user demodulator when 
there is no multiple-access interference. As can be seen 
from Fig 15 , the SER of the optimum multi-user demodulator 
is similar to the ideal value, that is, the lower limit of 
the SER of the single-user demodulator when there is no 
multiple-access interference . 

In Fig. 15 (and also in Fig. 16), the curve labeled 
"SINGLE-USER" was calculated in accordance with a 
theoretical expression descried in "Digital Communication" 
by J. G. Proakis (McGraw-Hill, 4th. ed., 2001, p. 261 
(hereafter, denoted as Ref. 18)). This value is similar to 
the practical lower limit in the range of E b /N 0 > 41og2 « 4.4 
dB. 

The performance of the multi-user demodulator varies 
depending on the relative value of the TH time t n k for 
respective users. As described above, the time hopping is 
performed to avoid collisions among users. If there is no 
pulse collision among users, the performance of the multi- 
user demodulator becomes similar to the ideal performance 
obtained by the single-user demodulator when there is no 
multiple-access interference. That is, in Fig. 15, the SER 
of the optimum multi-user demodulator indicated by the curve 
"ML" is similar to the lower limit of the SER of the single- 
user demodulator indicated by the curve "SINGLE-USER". 

On the other hand, when there is perfect collide 
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between pulses from different users, such as can occur when 
the TH times T x (1) and T/ 2> are both 0 and the pulse position 
of the user #1 and the pulse position of the user #2 become 
coincident with each other, the SER of the optimum multi- 
user demodulator represented by the curve "ML" in Fig. 15 
becomes identical to the SER of the single-user demodulator 
represented by the curve "CONV" . 

In the simulation shown in Fig. 16, it is assumed that 
the number of users, K, is equal to 4, and the symbol 
transmission interval T s is equal to 10 ns (and thus, the PRF 
becomes 100 MHz (= 1/10 (ns)). In the simulation shown in 
Fig. 16, the TH time t n k is given by a random number. 

In Fig. 16, an error floor occurs in the SER for both 
the optimum multi-user demodulator (indicated by a curve 
"ML") and the single-user demodulator (indicated by a curve 
"SINGLE-USER" ) . 

However, at E b /N G = 15 dB, the SER of the optimum multi- 
user demodulator (indicated by the curve "ML") is lower than 
that of the single-user demodulator (indicated by the curve 
"CONV") by a factor greater than 4. 

The above improvement in the SER by a factor greater 
than 4 is achieved because, in the multi-user demodulator, 
unlike the single-user demodulator in which demodulation is 
performed for each user, demodulation is performed jointly 
for all four users. The improvement in the SER achieved by 
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the multi-user demodulator , indicated in Fig. 16 , is very 
great . 

The processing sequence described above may be executed 
by hardware or software. When the processes are performed 
by software, a software program is installed on a general- 
purpose computer or the like. 

Fig. 17 illustrates an embodiment of the invention in 
which a program used to execute the processes described 
above is installed on a computer. 

The program may be stored, in advance, on a hard disk 
105 serving as a storage medium or in a ROM 103 which are 
disposed inside the computer. 

Alternatively, the program may be stored (recorded) 
temporarily or permanently on a removable storage medium 111 
such as a floppy disk, a CD-ROM (Compact Disc Read Only 
Memory), an MO (Magnetooptical ) disk, a DVD (Digital 
Versatile Disc), a magnetic disk, or a semiconductor memory. 
Such a removable recording medium 111 may be provided in the 
form of so-called package software. 

Instead of installing the program from the removable 
recording medium 111 onto the computer, the program may also 
be transferred to the computer from a download site via a 
digital broadcasting satellite by means of radio 
transmission or via a network such as an LAN (Local Area 
Network) or the Internet by means of wire communication. In 



this case, the computer receives, using a communication unit 
108 , the program transmitted in such a manner and installs 
the program on the hard disk 105 disposed in the computer. 

The computer includes a CPU (Central Processing Unit) 
102. When a user inputs a command by operating an input 
device 107 such as a keyboard or a mouse connected to the 
CPU 101 via a bus 101 , the command is transferred to the CPU 
102 via the input /output interface 110. In accordance with 
the command, the CPU 102 executes a program stored in the 
ROM (Read Only Memory) 103. Alternatively, the CPU 102 may 
execute a program loaded in a RAM (Random Access Memory) 104 
wherein the program may be loaded into the RAM 104 by 
transferring a program stored on the hard disk 105 into the 
RAM 104, or transferring a program which has been installed 
on the hard disk 105 after being received from a satellite 
or a network via the communication unit 108, or transferring 
a program which has been installed on the hard disk 105 
after being read from a removable recording medium 111 
loaded on a drive 109, thereby performing the process 
described above with reference to the flow charts and the 
block diagram. The CPU 102 outputs the result of the 
process, as required, to an output device such as an LCD 
(Liquid Crystal Display) or a loudspeaker via an 
input/output interface 110. The result of the process may 
also be transmitted via the communication unit 108 or may be 
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stored on the hard disk 105. 

In the present invention, the processing steps 
described in the program to be executed by a computer to 
perform various kinds of processing are not necessarily 
required to be executed in time sequence according to the 
order described in the flow chart. Instead, the processing 
steps may be performed in parallel or separately (by means 
of parallel processing or object processing). 

The program may be executed either by a single computer 
or by a plurality of computers in a distributed fashion. 
The program may be transferred to a computer at a remote 
location and may be executed thereby. 

Note that the UWB demodulator 13 can deal with the 
received signal r(t) regardless of whether there is 
multiple-access interference. 

As described above, the communication device, the 
communication method, and the program according to the 
present invention have great advantages. That is, the 
correlation is calculated between the received signal 
corresponding to the signals transmitted from a plurality of 
communication terminals by means of UWB (Ultra WideBand) 
communication technology and pulses at possible positions in 
the transmitted signals from the respective communication 
terminals, and the received signal is demodulated into 
original data issued by the communication terminals on the 
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basis of the calculated correlation taking into account 
interference among the transmitted signals from the 
respective communication terminals. Thus, a high- 
performance UWB communication demodulator is achieved. 



